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Abstract. Insects use flight muscles attached at the base of the wings 
to produce impressive wing flapping frequencies. The maximum power 
output of these flight muscles is insufficient to maintain such wing oscil¬ 
lations unless there is good elastic storage of energy in the insect flight 
system. Here, we explore the intrinsic self-oscillatory behavior of an 
insectile wing model, consisting of two rigid wings connected at their 
base by an elastic torsional spring. We study the wings behavior as a 
function of the total energy and spring stiffness. Three types of behav¬ 
ior are identified: end-over-end rotation, chaotic motion, and periodic 
flapping. Interestingly, the region of periodic flapping decreases as en¬ 
ergy increases but is favored as stiffness increases. These findings are 
consistent with the fact that insect wings and flight muscles are stiff. 
They further imply that, by adjusting their muscle stiffness to the de¬ 
sired energy level, insects can maintain periodic flapping mechanically 
for a range of operating conditions. 


1 Introduction 

The flapping of insect wings is a marvelous example of autonomous oscillations. In¬ 
sects use flight muscles attached at the base of the wings pQ to produce wing flapping 
frequencies that, in certain species, far exceed the animals’ neural capacity [2]. In these 
species, the contractile activity of flight muscles is maintained by a self-oscillatory 
mechanism that is under mechanical, not nervous, control [3J. Calculations of me¬ 
chanical power suggest that the maximum power output of flight muscles is adequate 
for the aerodynamic power requirements, but insufficient to also oscillate the wings’ 
mass unless there is good elastic storage of the inertial energy [3]. Elastic energy could 
be stored in several components of the insect flight system, including flight muscles. 
Insect flight muscles are very stiff such that, even at small operating strains, they can 
store elastically much of the inertial energy of the oscillating wing |3]. 

In this paper, we formulate and analyze the dynamics of an idealized insectile 
wing model, consisting of two rigid wings connected at their base by an elastic tor¬ 
sional spring (see Fig. [l]). Our goal is to examine the instrinsic self-oscillations of the 
wings under their own inertial load by drawing upon tools from nonlinear dynam¬ 
ical systems. This approach is aimed at complementing past and ongoing research 
efforts on flapping insect wings. The impressive aerodynamic performance of insects 
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Fig. 1 . Insectile wing model: two rigid wings connected by a torsional spring of stiffness k. 
The wings are free to flap (angle a) and rotate about their base point (angle 6). 


have stimulated a great deal of interest among biologists, physicists and engineers, 
whose focus is to decipher the biomechanics underlying insect flight and to translate 
this knowledge into design principles for engineered micro-air vehicles. Here, we give a 
brief overview of such work, with particular bias towards the aerodynamics of flapping 
flight. Early attempts at explaining the high-lift mechanisms in insect flight use quasi¬ 
steady aerodynamic analysis |4 . This approach provides valuable qualitative insights 
into the mechanisms of force generation in flapping flight but underestimates un¬ 
steady aerodynamic effects. The unsteady mechanisms responsible for the generation 
of remarkable lift forces in flapping flight has been examined in recent experimen¬ 
tal |5|6|7|8|9|10| and theoretical [11I12I13I14I15I16I17I18] studies, mostly emphasizing 
the importance of leading-edge and wake vorticity in force production |19l20j . Atten¬ 
tion turned recently to the stability of flapping flight in response to environmental 
disturbances [21|, with conflicting accounts of intrinsic instability [22123) and passive 
stability [24125126] . Given that an assessment of the passive stability of live organisms 
is not feasible experimentally, proxy models of inanimate flyers are proposed in [27128] 
and their passive stability is discussed in |29l30| . 

The present paper does not focus on wing aerodynamics, but on the intrinsic 
nonlinear dynamics of insect-like wings when subject to their own weight and elastic 
energy storage. Evolution and natural selection have made insect wings stiff enough 
to withstand the aerodynamic load m- Unlike birds and bats which have active 
muscles in their modified forelimbs to control the wing shape, insect wings twist and 
camber due to their elastic properties [32] . Wing flexural stiffness varies along the 
wing [33134] . However, we consider here a “lumped” spring model, where the wing 
and muscle stiffness are both modeled via an elastic torsional spring at the based 
of two rigid wings. We analyze the stability of this insectile wing system around 
the vertically down position and explore its global behavior as a function of total 
energy and spring stiffness. Interestingly, we identify regions of rotational motion 
(rolling end over end), chaotic behavior, and periodic flapping. The region of periodic 
flapping decreases as energy increases, and increases as spring stiffness increases. We 
conclude by discussing the potential implications of these findings on flapping flight. 


2 Model 

We consider two rigid wings, of mass m and length l each, connected at their base by 
a torsional spring of stiffness k (see Fig. [l]) . We let 2 a denote the opening angle of the 
wings, or the “flapping” degree of freedom, and 2a 0 denote the spring rest angle such 
that the internal spring torque is equal to 2k( a — a Q ). Further, the wings are free to 
rotate about their base point; their orientation angle from the vertical is denoted by 
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Fig. 2. Uncoupled wings (k = 0): Poincare sections at 0 — 0 for (a) E — E\ + E 2 = 2/3 and 
(b) E — 2. Insets show projections of the Poincare sections onto the plane 0 = 0. As energy 
increases, the compact energy surface becomes periodic in a. 


0. Ignoring aerodynamic forces, and under gravitational effects only, the equations of 
motion governing the behavior of a and 0 are given by 


(^ml 2 )e = - 
(| ml 2 )a = - 


mgl cos<a sin# 

mgl sin a cos 0 — 4^(<a — 


a 0 ). 


( 1 ) 


Here, g is the gravitational constant. In order to re-write the above equations in non- 
dimensional form, we use the length scale Z* = Z and time scale t* = y/2Z/3g, and 
introduce the dimensionless spring constant k* = 4 K/mgl. To this end, we get 


6 = — cos a sin 0 
a = — sin a cos 6 — k( a — a Q ), 


( 2 ) 


where we dropped the ()* notation with the understanding that all quantities are now 
non-dimensional. The total energy E of the wings is conserved, 

E = ~(<^ 2 + O 2 ) — cos a cos 6 + — a Q ) 2 . (3) 

For zero spring stiffness, k = 0, this wing model is equivalent to a system of two 
uncoupled rigid pendula. This system is integrable because it admits two conserved 
quantities, namely, E\ = + 0) 2 — cos(<a + 0)) and E 2 = — 6) 2 — cos(<a — 

#)), such that E = Ei + E 2 - Poincare section in the four-dimensional phase space 
(#, <a, a) is taken at 0 = 0. A depiction is shown in Fig. [2] for two total energy levels 
E = 2/3 and E = 2. At lower E , the energy surface is compact and closed, while at 
higher E, the energy surface is periodic in the (^-direction. 

For infinite spring stiffness, k = 00 , the two wings behave as one rigid pendulum, 
whose dynamics is governed by 0 = — cossin#. The frequency of oscillations of the 
pendulum depends on its shape, defined by the fixed parameter a. As a increases 
from 0 to 7 r, the oscillation frequency first decreases then increases. It reaches zero at 
a = 7t/ 2, where the gravitational torques on the two wings of the rigid pendulum are 
perfectly balanced and the system is in a state of equilibrium for all 0. 


3 Results 

We examine the behavior of the coupled insectile wings in Eq. © for finite spring 
stiffness k and rest angle a Q . Given the left-right symmetry of the wings, it suffices 
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Fig. 3. Insectile wings: (a) Stable flapping angle a* as a function of the spring coefficient k 
and rest angle a 0 . (b) Projection onto the (K,a 0 ) plane, showing the two stable regions for 
#* = 0 and #* = 7r, separated by a Q = 7t/2. 


to examine the dynamical behavior of a and # in the range [0,7r]. In this range, the 
equilibrium points (<a*,#*) of Eq. © can be identified as a function of k and a Q as 
follows 


* 

II 

o 

sum* = — a Q ), 


b 

II 

* 

sum* = K(a* — a Q ), 

(4) 

a* = 7t/2, 

cos#* = K,(a 0 — 7t/2). 



We study the linear stability of these equilibrium points by linearizing Eq. ([2]) about 


(— cos cos #* sum* sin#* \ /#\ , ^ 

y sin<a* sin#* — cos<a*cos 0* — k J ya J ' ' 

with the understanding that # and a in Eq. © are small. Linear stability analysis 
of the equilibrium points in Eq. 0 reveal the existence of two sets of linearly stable 
points, 



_ 7T 7T 1 

Fot-<a„<- + -, 

7T 1 7T 

F ° r 2’ 


9* = 0 and a* G [0, ^], 
9* = n and a* G [^> 7r ]- 


(6) 


Fig ,[3ja) illustrates the stable flapping angle a* as a function of k and a Q . Blue and red 
colors represent acute (a* < tt/2) and obtuse (a* > tt/2) flapping angles, respectively. 
A projection onto the (/c, a G ) parameter space is shown in Fig. [3Fb). Insets show the 
physical configurations of the wing system in these equilibrium positions. While the 
two configurations are the same, the internal torques used to balance the gravitational 
effects have opposite signs, namely, sum* = =F^(<a* — <a 0 ), where the minus sign 
corresponds to #* = 0, meaning that the internal spring is in a state of compression. 

An examination of the nonlinear response of the system around one of these 
equilibrium points shows that it exhibits regular or chaotic behavior depending on 
initial conditions (#,cq#,<a). Fig. hlJ shows examples of such behavior for E = 2, 
k = 4 and a Q = 7t/4, and two sets of initial conditions, namely, (7 t/36, 7t/3, 0.88, 2.5) 
and (7 t/36, 7t/3, 2.66, 0). To characterize the global nonlinear dynamics of the coupled 
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Fig. 4. Periodic (a) and chaotic (b) behavior around the equilibrium point (0, a*) shown as 
a red dot for E = 2, k — 4 and a Q — 7t/4. 
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Fig. 5. Poincare sections for k = 4, a Q = 7t/4 and E — 2/3, 2,10/3. 


wings, we construct Poincare sections (see, e.g., [35]) at 0 = 0, which reduces the 
phase space to a three-dimensional space (a, 0, d). Conservation of energy restricts 
the dynamical response of the system to lie on a two-dimensional energy surface 
within this space defined by ^(d 2 + 0 2 ) — cos<a + \k{ol — a Q ) 2 — E — 0. 

Fig .|5][a) depicts Poincare sections for k — 4 and three energy levels E = 2/3, 2,10/3. 
On each energy level, three regions can be distinguished: region I is characterized by 
periodic rotations; region II shows chaotic behavior; and, region III is characterized 
by periodic oscillations, which correspond to regular flapping behavior of the wings. 
Note that as energy increases, regions I and II grow, while region III shrinks. 

We now explore the effects of spring stiffness k on the nonlinear response of the 
wings. Fig. [6] depicts the energy surfaces for a constant energy value E — 2 and three 
different values of spring stiffness, k = 0.8,4 and 20. For small k, the behavior of 
the wings is predominantly chaotic. As k increases, the energy surface turns oblate 
in the ^-direction, indicating smaller flapping amplitudes, and the chaotic region II 
shrinks. A few comments on the limits as k goes to 0 and oo are in order here. At 
k = 0, the dynamical response is regular as depicted in Fig. [2j but for small non-zero 
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Fig. 6. Poincare sections for E = 2, a 0 = tv/ 4 and ft = 0.8,4, 20. The case of ft = 4 is same 
as the one shown in Fig. i(£ = 2). 


ft, the behavior becomes immediately irregular. That is to say, weak elastic coupling 
between the wings induces chaotic flapping. The limit as ft 0 is singular. Stronger 
elastic coupling between the wings tends to suppress the chaotic behavior, in favor 
of regular flapping oscillations which decrease as ft oc. In this non-singular limit, 
region I expands and shifts towards the rest angle, indicating small periodic flapping 
about the spring rest position a Q . At ft = oo, this periodic flapping is suppressed and 
the wings behave as a rigid pendulum, as mentioned in Section [2j 


4 Discussion 

We proposed an idealized two-dimensional model of insect-like wings. The model 
consists of two rigid wings connected at their base via an elastic torsional spring. We 
studied the passive dynamics of this system under gravitational effects, emphasizing 
the coupling between the wings flapping and rotational motions. We identified regimes 
where the wings flap stably and others where the wings rotate (roll end over end) 
or behave chaotically. For a given energy level, the flapping region increases as the 
stiffness of the torsional spring increases. These findings are qualitatively consistent 
with insect wings. Insect wings and the flight muscles at their base are known to 
be very stiff; even at the small operating strains, they can store elastically much of 
the inertial energy of the oscillating wings More importantly, our results suggest 
that, by manipulating the stiffness of their flight muscle, insects can maintain periodic 
flapping when operating at a range of energy levels. It is known that muscle force can 
be modulated using a number of mechanisms such as changing the muscle stiffness or 
length (introducing pre-tensioning in the muscle) [36] . We conjecture that modulation 
of muscle stiffness helps insect wings operate at a range of aerodynamic loads. 

The effect of aerodynamic forces on the wings dynamics can be accounted for using 
a vortex sheet model in the inviscid fluid context, as shown in Fig. [7] Here, the wings 
are modeled as a bound vortex sheet that satisfies zero normal flow. A point vortex is 
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Fig. 7. Elastic wings of spring stiffness n = 80 and rest angle a 0 = 7r/4 in rest fluid, 
(a) Snapshots of vortex shedding at four time instants highlighted in (b) and (c). (b) Time 
evolution of the orientation which approaches the vertically down position 6 = 0. (c) Flapping 
response a versus time. Initial conditions are set as 6 — 7r/6,a: = ir/3,0 — 0 and a = 0. 
Dissipation time Tdi ss = 1.22. 


released at each time step from the two outer edges, and the shed vorticity is modeled 
as a regularized free sheet |3Hsxl:ffll4()l41132| . No separation is allowed at the apex. 
We followed the algorithm in [38] for imposing the Kutta condition that determines 
the amount of circulation shed from the outer two edges at each time step. The vortex 
sheet model depends on the regularization parameter for the free sheet, which we set 
to S/l = 0.1. To emulate the effect of viscosity, we allowed the shed vortex sheet 
to decay gradually by dissipating each incremental point vortex after a finite time 
Tdi ss from the time it is shed in the fluid. This computational scheme is validated 
and used in 13QJ- Fi §- III shows a depiction of the wings’ flapping behavior, under 
both gravitational and aerodynamics forces and torques. Interestingly, the wings flap 
stably about the vertical position, even though the initial orientation 0 of the wings is 
perturbed away from the vertical position. Note that this stable flapping is observed 
at a larger value of spring stiffness than the ones explored in Fig. [6] A larger stiffness 
is needed to support the aerodynamic forces and torques due to the surrounding fluid. 
A detailed study of the aerodynamic effects on the wing behavior will be the subject 
of a future study. 

We conclude by noting that, in addition to neglecting aerodynamic effects, we 
made a number of assumptions in this paper to ensure the problem is tractable analyt¬ 
ically and computationally. We assumed the motion is planar, but three-dimensional 
wing rotations play an important role in force production in insect wings [516I2T]. 
We “lumped” all elastic components into a single torsional spring that couples the 
dynamics of the two wings. The elasticity of insect flight systems is distributed along 
the wings and in the flight and thoracic muscles Future extensions of this 

work will interrogate the effect of each of these simplifying assumptions on the wing 
flapping dynamics. 
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